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Growing and studying large Coulomb crystals, composed of tens to hundreds of thousands of 
ions, in linear quadrupole ion traps presents new challenges for trap implementation. We consider 
several trap designs, first comparing the total driven micromotion amplitude as a function of location 
within the trapping volume; total micromotion is an important point of comparison since it can limit 
crystal size by transfer of radiofrequency drive energy into thermal energy. We also compare the axial 
component of micromotion, which leads to first-order Doppler shifts along the preferred spectroscopy 
axis in precision measurements on large Coulomb crystals. Finally, we compare trapping potential 
anharmonicity, which can induce nonlinear resonance heating by shifting normal mode frequencies 
onto resonance as a crystal grows. We apply a non-deforming crystal approximation for simple 
calculation of these anharmonicity- induced shifts, allowing a straightforward estimation of when 
crystal growth can lead to excitation of different nonlinear heating resonances. In the anharmonicity 
point of comparison, we find significant differences between the trap designs, with an original rotated- 
endcap trap performing better than the conventional in-line endcap trap. 



I. INTRODUCTION 

Ion trapping by radiofrequency (rf) electric fields was 
first demonstrated in 1954 by Paul in a trap with hy- 
perbolic electrodes [1 . Attempts to maximize the field- 
free trapping volume led to the development of the linear 
Paul trap by Prestage in 1989 [2]. Although less har- 
monic than hyperbolic traps, the field- free central axis of 
linear traps often make theirs the preferred geometry for 
experiments involving multiple laser-cooled ions, includ- 
ing precision spectroscopy [3l [4], quantum information 
processing |^, and cavity quantum electrodynamics 
applications [71. 

When the thermal energy of a trapped ion cloud is 
cooled sufficiently below the repulsive Coulomb interac- 
tion energy, the ions settle into a Coulomb crystal with 
geometry defined by the trapping potential [8[ |9] . Crys- 
tallization of certain atomic ion species with closed-cycle 
electronic transitions can be accomplished directly by 
laser cooling [lQl[TT]; crystallization of other species (e.g., 
molecular ions) can be accomplished sympathetically by 
laser cooling of a co-trapped atomic species [12 or in 
situ formation from a pre-cooled reactant [13]. Crystal- 
lization can be prevented by various heating mechanisms. 
Two such heating mechanisms, micromotion heating and 
nonlinear resonance heating^ become more problematic 
as a crystal grows, due to sampling of larger rf fields 
and of nonharmonic regions of the trapping potential, 
respectively. Thus, these heating mechanisms can limit 
the minimum temperature obtained for a given crystal 
size and also the maximum obtainable crystal size. Our 
work is motivated in part by recent difficulties growing 
large Coulomb crystals in a non-standard trap geometry, 
leading to speculation that axial rf fields near the endcaps 
(fringing fields) might limit crystal size [T4|. 

In this paper, we compare the suitability of four linear 
Paul trap designs for growth and study of large Coulomb 
crystals. First, we compare the total micromotion (the 



unavoidable driven motion at the rf frequency) which 
results in non-resonant transfer of rf energy into ther- 
mal energy of the trapped sample (micromotion heat- 
ing). Second, we compare the component of micromo- 
tion along the trap axis (the axis of smallest micromo- 
tion) , which induces first-order Doppler shifts in precision 
spectroscopy on large Coulomb crystals. Finally, we com- 
pare trap anharmonicity, which can lead to excitation of 
nonlinear heating resonances as a crystal grows and its 
normal mode frequencies are shifted onto a resonance. 
To characterize these anharmonicity-induced frequency 
shifts, we use a simple non-deforming crystal approxima- 
tion, which allows a straightforward estimation of when 
crystal growth can lead to excitation of different nonlin- 
ear heating resonances. 



II. TRAP DESIGNS 

A linear quadrupole ion trap (Fig. [T]) is formed from 
four parallel cylindrical electrodes of radius re, held at 
a separation ro, with an rf voltage signal applied to the 
electrodes. The resulting electric field confines the ion 
motion radially; axial confinement (in the 2;-direction) is 
provided by endcap electrodes (not shown in Fig. [T]), to 
which a static voltage is applied. 

Near the trap center, the electric potential is given by 

^0 



where Vrf and Vtrf are the voltage and frequency of the 
applied rf drive, Vec is the static voltage applied to the 
endcap electrodes, zq is half the distance between the 
endcap electrodes, and is a geometric factor. While this 
description is accurate near the trap center and along the 
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FIG. 1. End view of a linear quadrupole trap. A sinusoidally 
varying voltage is applied to the rod electrodes, with one pair 
held 180° out of phase with the other. 



z-axis, significant deviations can occur elsewhere in the 
trap interior. 

The radial motion of an ion in the potential of Eq. ([T]) 
is described by the Mathieu differential equations, and 
the stability of this motion is expressed using Mathieu 
parameters which depend only on Vrf^ Vec^ ^r/, and 
the charge-to-mass ratio q/m. Solutions to the Mathieu 
equations contain regions of parameter space where ion 
motion is stable [15], and trap parameters are chosen to 
operate in one such stable region. 

The motion of a trapped ion can be described approx- 
imately as a superposition of micromotion, in which the 
ion oscillates at the rf frequency, with secular motion, in 
which the ion oscillates in a time-independent pseudopo- 
tential (j)rf{x^y^z) at slower secular frequencies ujx^ 
and ujz' An analytic relationship exists between the time- 
dependent rf potential and the time-independent pseu- 
dopotential [16] : 
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The endcap electrodes generate a potential ^ec, and the 
total potential governing secular motion is 



l^trap 



Prf 



(3) 



We compare four different linear trap implementations 
(Fig. [2j with corresponding voltage given in Table [l|): 
a conventional in-line endcap design with two different 
voltage configurations (designs A and 5), a plate end- 
cap design (design C, similar to a design analyzed in 
Ref. [17 ), and an original rotated endcap design (design 
D). In all trap designs, zq = 8.5 mm, ro = 4 mm, and 
Ve = 4.5 mm, with ro and rg chosen to agree with the 
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FIG. 2. Trap designs analyzed in this paper. Designs A and 
B share the same geometry, but differ in applied voltages. In 
the offset view at right, designs C and D are shown with the 
nearest plate endcap pulled away from the rest of the trap to 
allow the trap interior to be seen. 



TABLE I. Operating voltages applied in different trap de- 
signs. Here V(t) = -\-VrfCos(Qrft). 
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Vz 




A 


V(t) 


-v{t) 






B 


V{t) 


-V{t) 


Vea + V{t) 


Vec - Vit) 


C 


V{t) 


-V(i) 


Vec 


N/A 


D 


V{t) 


-V(i) 


Vic 


N/A 



optimal ratio re=1.14511ro for minimizing the leading- 
order contribution to anharmonicity in the rf potential 
[18] . Design C uses thin circular plates of radius 2re + ro 
as endcap electrodes, with a 2 mm radius hole left for 
axial optical access. Design D has four small cylindrical 
endcap electrodes of radius 1.75 mm, rotated 45° around 
the z-axis; the central axes of these endcap electrodes 
are on the same radius as the central axes of the rf elec- 
trodes. Our group has built a functional trap whose key 
geometric features are shared with design currently 
being used to trap ^^^Ba+ ions. All numerical simula- 
tions in this paper use values corresponding to ^^^Ba+ 
{q = 1.602 X 10-1^ C, m = 2.292 x 10"^^ kg) and our 
operating rf frequency {Vtrf = 27r x 3.00 MHz). 



III. COMPUTATION OF TRAPPING 
POTENTIALS 

Our metrics of trap design comparison, presented in 
Sec. IV and Sec. |VI[ are calculated from the trap po- 
tentials and their gradients. The potentials (j)ec and (prf 
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TABLE II. Vrf and Vec used in comparisons and calculated 
value of K for each trap. 



Trap Design 


Vrf (volts) 


Vec (volts) 


K 


A 


800 


5.0 


0.15 


B 


800 


5.0 


0.15 


C 


800 


2.7 


0.27 


D 


800 


51 


0.014 



are found using the finite element method to solve the 
Laplace equation, with the trap geometries discretized 
by the meshing software Gmsh [19]. 

To allow evaluation of the potentials at locations not 
on the mesh, the numerical solutions are fitted to an ex- 
pansion in associated Legendre polynomials 

Imax Imax 

{cos(m(b) for m > 
^ ^ ~ (4) 

sin{m(j)) for m < 0. 

We exclude terms with odd values of / and m from our 
fitting due to symmetry. 

An analytic form of the gradient of Eq. Q is straight- 
forward to find using 

dPl^{cos 6) _ 

lcosOP^{cos 6) - {l^m)P%cos 0)) 
Vl - cos^O ' 



(5) 



which can be derived from a more general identity given 
by [20]. Fitting the Laplace solutions to the above ex- 
pansion thus permits efficient evaluation of (pfu or its 
gradient V(l)fit at arbitrary location from the coeffi- 
cients {Cim}' In this way, we determine (j)rf and (pec 
up to the scaling set by choosing Vrf and Vec Our 
choices of these voltages, shown in Table [ll] are made 
to equalize the single particle secular frequencies ujz and 
ujr across all designs at experimentally reasonable values 
ujx = uoy = 271 X 418 kHz and cj^ = 27r x 18.9 kHz. 

Single particle secular frequencies are calculated by 
evaluation of (j)rf and 0ec along a radial or axial trace 
through the origin from r=— ^ to r=+^ or z=—^ to 
z=^ respectively. The resulting values fit a polynomial 
whose quadratic coefficient rapidly converges with finer 
discretization of the domain or increase of the polynomial 
fitting order. From this coefficient C2, we also find the 
values oi K = ^J2qC2/m in Table |ll| 

The effect of inaccuracy in the Laplace solutions is esti- 
mated by observing convergence of calculated values such 
as K with increasing meshing density, which is found to 
cause variation in all reported results of less than a few 
percent unless otherwise noted as we approach the high- 
est mesh densities eventually used - of order 10^ nodes 
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FIG. 3. (Color online.) Contour plots of total micromotion 
amplitude, labeled in microns. Designs B and D are indistin- 
guishable at this scale. 



for all designs. We discuss the implications of the result- 



ing uncertainty in Sec. VII Tuning of other parameters, 
such as the fitting order in Eq. Q, is found to cause 
negligible variance by comparison. 



IV. MICROMOTION 

Although radial and axial micromotion vanish, respec- 
tively, along the trap axis and at z = 0, a large crystal 
necessarily extends beyond these regions. Micromotion 
heating results from the transfer of micromotion energy 
into secular energy, which leads to an elevated secular 
temperature. Although the scaling is non-trivial, simu- 
lations show that the micromotion heating rate strongly 
increases with micromotion amplitude and with secular 
temperature [2T1[22]. As a crystal grows, additional ions 
are held at locations with increasingly strong rf fields, 
causing them to experience increasingly large micromo- 
tion. The resulting heating limits the minimum obtain- 
able temperature for a crystal of a given size, and even- 
tually prevents further crystal growth. 

The electric field created by the rf electrodes may be 
written as 



^(x, y, z, t) = ^o(^, z)cos{Qrft)' 



(6) 



For an ion held in place in a stationary crystal, in the 
limit of small micromotion amplitude Eq may be taken 
as constant as the ion moves, and the micromotion of the 
ion is harmonic with amplitude 



A{x,y,z) 



\V(t)rf{x,y,z,t = 0)1 



(7) 



With Vrf as set in Table [H] we find (prfi^^y^z^t = 0) 
and calculate the resulting total micromotion amplitude 
A{x^y^z). The resulting micromotion contour is plotted 
m Fig [3| All desi gns are comparable in the amplitude 
of total micromotion for any realistic trapping volume. 



Interpretation of this result is discussed in Sec. VII 
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TABLE III. RMS average axial micromotion amplitudes over 
the cylindrical volume r < 500/im, \z\ < zo = 8.5mm. Listed 
values for A and C change by no more than several percent 
when Laplace solutions are calculated over an order of mag- 
nitude in grid density. Values for B and D change by up 
to tens of percent. For averaging, the volume is discretized 
by a Cartesian grid with h = 20/im, which is verified to be 
sufficiently dense. 



Trap Design 


RMS Axial Micromotion Amplitude (/xm) 


A 


1.5 


B 


0.04 


C 


2.9 


D 


0.08 



We also compare the size of the z component of mi- 
cromotion in each trap, which can limit the precision of 
spectroscopy. Due to the smallness of these values (sev- 
eral orders of magnitude smaller than the total ampli- 
tude) numerical convergence comparable to Fig. [s] is not 
observed. Rather than plot contours, we compute the 
root mean square (RMS) average over a large trapping 
volume for each design in Table |III[ producing a single 



quantity whose convergence is verified. See Sec. |VII| for 
discussion. 



V. NONLINEAR RESONANCE HEATING 

Nonlinear resonance heating can occur if a resonance 
condition is satisfied between secular frequencies and the 
rf drive frequency. For linear Paul traps, the condition is 



rixUJx + UyUJy + UzUJz 



rf 



(8) 



where n^^, n^y, and are integers [23 . The correspond- 
ing condition for a single particle in a hyperbolic trap has 
been derived [24], with the resonances weakening with 
larger n. This behavior has been confirmed in both hy- 
perbolic traps [25^ and linear traps [23 . For sufficiently 
cold trapped samples, it is understood that the frequency 
of the center of mass (COM) modes, or other normal 
modes, should be used in Eq. ([8| rather than the single- 
particle frequencies [26] . For an infinitely long linear Paul 
trap, or for a small crystal in a finite trap, where the ax- 
ial rf drive vanishes, a non-deforming crystal is only ex- 
pected to be heated on resonances with = 0. However, 
a crystal which is large enough to sample axial fringing 
fields near the endcaps is expected to be excited also by 
resonances with > 0. 

Anharmonicity in the trapping potential causes the 
mode frequencies to shift as a crystal grows and samples 
less harmonic regions of the trap potential. Eventually, 
the anharmonicity-shifted mode frequencies of a growing 
crystal will meet the resonance condition of Eq. (|8|, and 
some level of heating will occur, potentially halting fur- 
ther crystal growth. Even for a crystal of definite size. 



it is non-trivial to predict the precise response when it 
is swept through a given nonlinear heating resonance. 
Predicting the heating response in the non-equilibrium 
scenario of a growing crystal, as anharmonic frequency 
shifts cause it to cross a resonance, is yet more compli- 
cated. Modeling of the heating rates and whether cross- 
ing a given subharmonic resonance will actually melt a 
crystal or prevent further growth is beyond the scope of 
this analysis. 

Since it is difficult to determine whether crossing a 
heating resonance will limit crystal growth or cause other 
deleterious effects, in the current analysis, we take a 
conservative approach of attempting to design a trap in 
which a growing crystal will avoid low-order resonances 
altogether. We consider the anharmonicity-induced 
shifts of the COM frequencies for a non-deforming crys- 
tal; the non-deforming approximation is exact in the 
zero-temperature limit. In light of the resonant condi- 
tion of Eq. (|8|, considering the COM frequency shifts 
provides a simple estimate of when trap anharmonicity 
could limit crystal growth. Other modes of higher order 
than the COM mode can also be excited by nonlinear 
resonance heating, but the COM shifts are easily calcu- 
lable and set the scale for anharmonicity-induced shifts 
in higher modes. 

A previous analysis by the Marseille group (Pedregosa 
et. al. [12) specified a single quantity describing the 
radial trap anharmonicity averaged over a trapped sam- 
ple. The Marseille approach has the utility of providing a 
single figure of merit by which to compare different trap 
designs; however, it does not consider axial anharmonic- 
ity (potentially important for long crystals) or provide 
an estimate of what level of anharmonicity is acceptable 
in order to avoid heating resonances. The analysis we 
present here is similarly simple to calculate, and at the 
cost of not providing a single figure of merit, seeks to 
include axial anharmonicity and to provide a mechanism 
for determining whether a trap is sufficiently harmonic 
to safely work with crystals of some size. 



VI. COM FREQUENCY 

For a zero-temperature crystal in a perfectly harmonic 
trapping potential, the COM oscillation frequency is the 
same as the secular frequency of a single particle. How- 
ever, if the temperature of the crystal is nonzero or if 
the trapping potential is anharmonic, the crystal deforms 
during oscillation, and the COM frequency and single 
particle frequency are no longer equal. (Space charge 
shifts, which are well-understood in the limit of each ion 
moving independently in a background potential describ- 
ing the distributed charge of the other ions, occur at 
high temperatures but not for a non-deforming crystal.) 
Molecular dynamics (MD) simulations can find the COM 
frequency under finite temperature conditions, but this 
approach is computationally intensive for a large crys- 
tal. To understand the effects of anharmonicity with less 
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computational overhead, we introduce the approximation 
of a non-deforming crystal in an anharmonic trap. This 
approach is computationally simple enough to be useful 
in designing traps for large Coulomb crystal experiments. 

A single ion in a zero temperature crystal is held in 
place by the cancellation of two forces: that due to the 
trap potential Ftrap and that due to the Coulomb repul- 
sion of all other ions in the trap Fions • In the approxima- 
tion that the crystal does not deform during oscillation, 



remains constant while Fo. 



Ft 



trap 



is, to lead- 



ing order, a restoring force resulting in simple harmonic 
motion. 

Ftrap can be calculated from (ptrapi evaluated on a grid 
{xi,yj,Zk) where Xi = ih.yj = jh, and Zk = kh; i, j, k 
are integers; /i=50/im is the step size used in our calcu- 
lation. To illustrate our approach, we discuss calculation 
of the restoring force for small displacements in the z 
direction, at fixed Xi and yj. By geometrical symmetry 
^trap is an even function of z^ so the numerically calcu- 
lated values can be fitted to a polynomial series of even 
terms: 



(t^trap{^) 



Co + C2^2 



Xi,yj 



(9) 



The axial force on an ion in a potential of this form is 



trap 



(z) = -q{2C2Z + ... + lOCioz^). 



(10) 



We now consider an ion whose equilibrium location in a 
crystal is z = Zk. An expansion of Ftrap about Zk gives 

Ftrap{z) = Ftrap{Zk) + - Zk)Fi^^p{Zk) + ... (H) 

Since z = z^ is an equilibrium position for this ion, the 
Coulomb force due to all other ions Fions must exactly 
cancel the first term of this expansion. In the approxi- 
mation of a non-deforming crystal, Fions is a constant, 
and the total force on the ion is 



Ftrap H~ Fions 



= -q{2C2^...^90Ciozl){z-Zk) (12) 
kz (z-Zk), (13) 



where kz\xi,yj,zk is a local spring constant describing the 
strength of the axial restoring force on an ion about its 
equilibrium location {xi^yj^ z^). As seen from Eqs. (10) 
and (12), when 7^ 0, anharmonic terms of the trap 
potential contribute to the spring constant /c^^ . 

The restoring force on a crystal is the sum of the restor- 
ing forces on the individual ions. For a non-deforming 
crystal, the result is oscillatory motion at a frequency set 
by the averaged local spring constants, with the average 
taken over the ion locations. In our calculation we aver- 
age instead over the grid locations, for points interior to 
the crystal volume. In this way, we can estimate the axial 
COM frequency for a crystal, given (j)trap and a specified 
crystal geometry. We compute the COM frequencies for 
radial motion by the same method, again using a fitting 
polynomial of the same order. 
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FIG. 4. (Color online.) COM axial secular frequency 
as a function of ion number. Curves for Designs A and B 
are overlapping. Crystals are grown until their axial extent 
reaches ±zo/2; for equal length, design D grows an apprecia- 
bly smaller crystal. Ripples result from finite evaluation grid 
size. (See Sec. VII for discussion of both ion number in D as 
well as finite grid effects.) 



Fig. |4] shows the change in axial COM frequency in 
each trap design as crystal volume grows. Crystal extent 
is determined by evaluation of (j)trap 011 a Cartesian grid 
with spacing /i=50/im, with all points below a fixed cutoff 
included, and ion number is then determined by the ion 
density [27^ . 



^^^^ 



(14) 



where ujz and ujr = oOx = ojy are the single particle secu- 
lar frequencies. We take this density as a constant; the 
anharmonicity-induced position dependence is negligible. 
Implications of anharmonicity-induced frequency shifts 
for nonlinear resonance heating are discussed in Sec. |VII| 
Radial COM shifts are also computed, with fittings 
made to radial rather than axial traces of (ptrap- Over 
the volume occupied by the largest crystals in Fig. |4j no 
difference among the designs can be determined at our 
level of computational accuracy. By varying all calcula- 
tion parameters, we bound radial COM shifts at no more 
than tenths of a percent. 



VII. DISCUSSION 

No difference in micromotion heating rate is expected 
among the different designs, since total micromotion am- 
plitudes are comparable. Indeed, in all cases axial micro- 
motion amplitudes are much smaller than total ampli- 
tudes, indicating that rf-fringing should not limit crystal 
growth for these designs. The actual limit on crystal 
size due to micromotion heating is difficult to estimate, 
but it can in principle be found from molecular dynamics 
simulations. 

Axial micromotion, quantified by RMS averaging over 
a volume bounding any realistic crystal, is at least than 
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an order of magnitude smaller in designs B and D (in 
which the rf electrodes extend to z > zq) than in A and 
C (in which the rf electrodes end at z = zq). Since axial 
micromotion results from fringing of the rf fields, these 
differences are expected, given the different rf boundary 
conditions at z = ±zo. Given a conservative interpreta- 
tion of the observed convergence of this average, B and 
D may be said to have comparable axial micromotion, 
while the situation in A is much worse and in C worse 
stih. 

The size of COM frequency shifts during crystal 
growth, found here to differ among the considered de- 
signs, is an indicator for the effect of nonlinear resonance 
heating on trap operation. As the experimenter grows 
a crystal, adding ions continuously from an ion source, 
small COM frequency shifts (say a few percent) are not 
expected to present significant challenge. The COM fre- 
quencies are unlikely to shift onto a nonlinear resonance 
during the intermediate stages of crystal growth, and 
the experimenter can load the largest crystal possible 
before checking if different voltages increase the maxi- 
mum crystal size. For larger size-dependent shifts, trap- 
ping voltages could, in principle, be dynamically tuned 
to avoid resonances as the crystal grows, but dynamic 
tuning would become increasingly challenging for larger 
shifts. Thus the designs exhibiting large COM frequency 
shifts are expected to be more problematic for growing 
large crystals. 

As an example of these concerns, we note that for 
the frequencies considered in this study, excitation of the 
nonlinear resonance + |n^| = 7, = 3 is expected 
when the axial COM frequency increases by 5-10%, a 
shift which occurs in design D only for a crystal contain- 
ing tens of thousands more ions than in other designs. 
(Resonances of this order and higher have been observed 
experimentally in a linear trap [23 .) The expected abil- 
ity of D to load more ions than the other designs before 
exciting this resonance is a direct result of the smaller ax- 
ial COM shifts in D. In general, since the rotated endcap 
design D has the smallest axial COM shifts it is expected 
to be the least susceptible to nonlinear resonance heating. 

The plots in Fig. [4] terminate at crystals of equal ax- 
ial length. For a perfectly harmonic trapping potential, 
the aspect ratio of a crystal will be independent of size 
and equal to the ratio of the secular frequencies, which in 
these traps were tuned to be equal in the harmonic trap- 
ping region. Accordingly, eventual differences in aspect 
ratio are due to different degrees of anharmonicity sam- 
pled by the trapped ions. In Fig. [5j we plot the expected 
change in aspect ratio during crystal growth in each de- 
sign, calculated by a simple comparison of the heights of 
the radial and axial potentials along x and z respectively. 
The tendency of design D to maintain larger aspect ra- 
tio than other designs explains the ion numbers seen in 
Fig. [4j at equal length, the crystals in designs A, B and 
C are radially wider than D and so contain more ions. 
Intuitively, the tendency of D toward radially narrower 
crystals helps explain the smaller axial COM shifts in 
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FIG. 5. (Color online.) Aspect ratio as a function of crystal 
size, with aspect ratio estimated by comparison of the heights 
of the radial and axial potentials along x and z respectively. 
A is overlapped by B and not visible. 

since the regions nearer the z axis are expected to be 
more harmonic than those radially further away. 

The ripples in Fig. |4] caused by the finite size of the 
evaluation grid do not interfere with the conclusion that 
the COM frequency shift in D is smaller than in the other 
designs, as the trend is clear at much sparser grids, where 
the ripples are even more pronounced. The Laplace so- 
lutions are found to converge sufficiently to determine 
the COM shifts to within a few percent, and all other 
parameters-the associated Legendre and polynomial se- 
ries fitting orders and the evaluation grid density-are ad- 
justed until the resulting convergence in COM frequency 
shifts is not limiting. 



VIII. CONCLUSION 

We have compared the suitability of four trap designs 
for experiments on large Coulomb crystals. The com- 
parison is based on the desire to minimize micromotion 
heating and anharmonicity-induced nonlinear resonance 
heating. 

We have compared total micromotion as a function of 
position within the trap volume. Micromotion heating, 
related to the total micromotion amplitude, is expected 
to be comparable in all designs. 

We quantify axial micromotion by RMS averaging over 
a large trapping volume, finding that the resulting limit 
on spectroscopic accuracy in B and D to be much bet- 
ter than in either A or C, with C expected to be least 
suitable of all for precision spectroscopy on large crystals. 

In order to compare the susceptibility of the differ- 
ent trap designs to nonlinear resonance heating, we have 
used a non-deforming crystal approximation to map the 
anharmonicity-induced axial and radial COM frequency 
shifts as a function of ion number. Shifts of COM fre- 
quencies, as well as other normal mode frequencies for 
which COM shifts set the scale, are of concern since a 
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growing crystal in an anharmonic trap will eventually 
cross a nonlinear heating resonance, leading to instan- 
taneously elevated temperatures and potentially halting 
further crystal growth. We find that a novel rotated end- 
cap trap has the smallest axial COM shifts and observe 
that the relative stability of aspect ratio with increasing 
crystal size in this trap makes the result intuitive. We 
suggest that traps aiming to grow large Coulomb crystals 
ought to minimize axial COM shifts, and we provide a 
method for quantifying these shifts. 

The analysis presented here is simple to implement and 
should be predictive of trap behavior, so it is especially 
useful for taking a conservative approach to designing a 
trap for experiments on crystals of a desired size. How- 
ever, our analysis is not capable of predicting precisely 
when crystal growth will be limited by micromotion or 



by a given nonlinear heating resonance being crossed dur- 
ing crystal growth. To make detailed predictions of these 
limitations on ultimate crystal size and temperature, MD 
simulations could be used. It would also be quite inter- 
esting to experimentally investigate the upper limits on 
crystal size set by anharmonicity-induced nonlinear heat- 
ing resonances. 
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